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THE  AUTOCORRELATION  AMD  JOINT  DISTRIBUTION  FUNCTIONS 

OF  THE  SEQUENCES,  {-  j2),  (^(j+T)2} 

_ _  m  m 


by 

D.  L.  Jagerman 


1.  INTRODUCTION 

The  present-day  extensive  use  o£  Monte  Carlo  procedures  necessitates  the 
careful  Investigation  of  methods  for  the  generation  of  random  numbers.  In  its 
simplest  form,  the  underlying  principle  of  all  Monte  Carlo  procedures  finds  its 
expression  in  the  following  theorem. 1 

THEOREM:  Let  (x^)"  he  «  sequence  equidistributed  over  (0,1),  and  let  f(x)  be  a 
function  Riemann  integrable  on  (0,1);  then 
1 

N  n 

2  f(x  )  ~  Nj  f(x)dx  . 


Thus,  the  theorem  states  that  sample  averages  approximate  the  value  of  an  in¬ 
tegral.  It  is  to  be  noted  that  the  only  property  of  the  sequence  (x^)”  employed 
is  its  equidistribution. 

In  applications,  one  may  employ  several  equidistributed  sequences  simul¬ 
taneously;  accordingly,  new  requirements  may  arise.  The  sequences  may  be  employed, 
for  example,  as  bases  for  decision,  in  which  case  it  may  be  required  that  they 
be  independent.  ThUs;  depending  on  the  Monte  Carlo  problem  considered,  equi¬ 
distributed  sequences  may  be  required  to  possess  other  random-number  character- 
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istics.  Let  the  sequences  (x^)*,  <x^+T>“  be  designated  reapectlvely  by  x, 

in  which  t  is  a  non*zero  integer,  then  an  additional  desirable  character* 
istic  is  the  statistical  independence  of  x,  A  sequence  exhibiting  such 

characteristics  is  ((aj  })^  in  which  a  is  an  irrational  number.  The  symbol 
(x)  is  employed  to  designate  the  fractional  part  of  x. 

However,  from  the  viewpoint  of  the  practical  utilization  of  the  sequence 
suggested  above  by  means  of  a  digital  computer,  it  is  necessary  to  replace  a 
by  a  rational  number,  and,  hence,  to  lose  some  of  the  precision  with  which  the 
characteristics  discussed  above  are  satisfied.  Accordingly,  the  sequences  which 
will  be  studied  are 

*  -  «*  j2c1-  «(,)  ■  «>o2»ri  • 


The  integers  a,  m  are  taken  relatively  prime.  Of  particular  interest  will  be 
the  deviation  of  the  characteristics  of  these  sequences  from  the  ideal  random* 
number  characteristics. 

Let  p(x)  be  given  by  p(x)  *  -  (x) ;  then  the  autocorrelation  function 

00 

KT)  of  a  sequence  (Xj)^  Is  defined  by 

*C0  -  £  jiiPfrj)  P<xJ+T)  • 

For  the  sequence  to  be  studied,  this  takes  the  form 

-  i  ”£Sp<5  ft  • 


I 
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It  will  be  shown  that,  uniformly  in  1  for  the  range 
1  <  t  <  v/m  , 


the  autocorrelation  function  is  small;  that  is,  quantatively, 

_  1 

U(t)|  <  m  2[|(2  +sT2)v(m>  In2  nH-30  •  2V*m)  In  m]  , 


provided  (a,m)  *  1  and  m  >  36.  The  function  v(m)  is  the  number  of  distinct 

s  2  ni*  X 

prime  divisors  of  m.  Thus,  the  sequence  ({—  j*})  is  approximately  uncorrelated. 
Let  Hq(x)  be  given  by  Hq(x)  ■  a  +  p(x)  -  p(x-a);  then  the  joint  distribution 

(1) 

function  G(a,0)  of  the  sequences  x,  xv  '  is  given  by 

5  V’V  • 


8  2  m-1 

For  the  sequence  ({—  j*})  ,  this  takes  the  form 


:  ISjvi  i2>  V:<J«>2> 


It  will  be  shown  that,  uniformly  in  t  for  the  same  range  stated  above,  and  under 
the  same  conditions  on  a,  m, 

_  1 

|  G(a,{3)  -op |  <  m  2  (2  +V*2)v(m)  (|  In2  mf88  In  m)  . 

Thus,  the  sequences  ({—  j2})™  *,  ((~(j+i) 2))™*  *  are  approximately  independently 
equidistributed  over  (0,1). 
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The  above  enumerated  properties  show  the  possible  applicability  of  the 

sequences  ((—  as  random  numbers  in  Monte  Carlo  procedures, 

m  o  m  o 

However,  an  important  question  is  the  behavior,  from  the  viewpoint  of  random* 
number  characteristics,  of  consecutive  portions  of  the  complete  sequences.  This 
is  being  studied  by  the  author,  and  an  investigation  of  the  question  will  appear 
in  another  paper. 


2.  ANALYTICAL  DISCUSSION 

The  proofs  of  the  main  theorems  require  the  establishment  of  several  lemmas. 
Lemma  1;  t  >  1,  |q|  <  1 


— >> 


p00 


-  L 

1  <  h  <  t 


sin  2nhx 
nh 


+  q  min  (1, 


) 


in  which  ||x||  denotes  the  distance  from  x  to  the  nearest  integer.  It  is  under¬ 
stood  that  when  x  is  an  integer,  the  estimate  1  is  used. 


Proof:  The  Fourier  series  for  p(x)  is 


P(x) 


«  sin  2nhx 
h-1  nh  ' 


(1) 


hence,  it  may  be  necessary  to  establish 


£ 

h>t 


sin  27thx 
nh 


<  min  (1, 


2«t  x 


(2) 


The  following  standard  theorem  derived  from  Abel's  transformation  of  series  will 


be  used: 
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'll  ^  o, 


£ 

pSMbp 


<  B 


[•M* ib£  -  *M  8 


(3) 


Also  standard  is  the  following  estimate: 


sin  2xpx 


(4) 


Applying  the  inequalities  of  Equations  (3)  and  (4),  one  has 


£ 

h>t 


sin  2-rhx 
xh 


-  2rr([t]+l)l|x 


< 


1 

2nt||x||  ‘ 


(5) 


If  2xt||x||  >  1,  then  Equation  (2)  has  been  established.  Consider  now  the  case 
2xt||x||  <  1;  then 


£ 

h>t 


sin  2nhx 
xh 


?  sin  2xhx  _  „  sin  2xhx 

h-1  xh  1  <  h  <  t  flh 


(6) 


Thus 


sin  2xhx 

<  1  . 

sin  2xhx 

L 

h>t 

xh 

-  2  + 

1  <  h  < 

t  xh 

<1  + 

£ 

Isin  2fthllxlil 

-  2  • 

L  <  h  <  t 

xh 

(7) 

+  2||*||t  <i  +  i< 


1 


Equation  (2)  is  now  established.  The  Fourier  series  for  p(x)  does  not  equal 
p(x)  when  x  is  an  integer;  however,  with  the  understanding  stated  in  the  lemma, 
the  lemma  remains  correct  also  in  this  case. 
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Lemma  2: 


2  sin  2tthx 
1  <  h  <  t  «h 


<! 


Proof: 


One  has 


sin  2ithx  _  “  sin  2nhx  _  z  sin  2nhx 
1  <  h  <  t  *h  h»l  *h  K>t  *h 


(8) 


Hence, 


2  sin  2-hx 
1  <  h  <  t  «h 


*i* 


_  sin  2*hx 
t£t  «h 


<1 


(9) 


The  inequality  of  Lemma  1  was  used. 


Lemma  3:  t  >  1,  U|  <  1 


P(x)p(y)  -  E  E  sin  2ithx  •  sin  2n^r 

PWPW  1  <  h  <  t  1  <  1  <  t  n2hi 


1 TlUin  (l>  +  min  <l'  ^R)] 


Proof:  Use  of  Lemma  1  yields 


P(x)p(y) 


E 

1  <  h  <  t 


E 

1  <  l  <  t 


sin  2nhx  •  sin  2n/y 
«2h  l 


+ 


n 


E 

1  <  h  <  t 


sin  2rhx 
jth 


min  (1, 


1 

2*t||y|| 


■)  + 


E 

1  <  i  <  t 


8 in  2ttiy 

nl 


(10) 


min  (1, 
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2  1  1 
*1  min  min  ds  ?  • 


2*t  x 


2*  1 1|  y  || 


Lemma  2  allows  one  to  write 


p(x)p(y)  -  £  2  al-S-M*.-  ,aln  2«iy  + 

1  <  h  <  t  1  <  l  <  t  *2hi 


I  ^min  <l>  jdbr)  +  min  C1,  ^rrnr)l  + 


2«t||y||J 


(U) 


T)  min  (1,  “  ■  || 1  ,,)  min  (1,  „_^n_j|) 


2xt  x 


2*t||y|| 


Observing  that 


m*n  (1,  „__ii__n)  ®ln  (1,  ^_^u„i|)  £  mln  (1,  2nt||y^  — 


2*t||> 


2«t||y|| 


“i11  (i*  o_ *. M __ii )  +  min  (i,  r-Tirjr) 


2*t  x 


2«t||y|| 


(12) 


the  lemma  follows. 

l^t  f(j),  g(J)  be  given  functions  of  the  integral  variable  j.  Define  e(x) 
by 


,  s  i2nx 
e(x)  -  e  , 


(13) 


R  by 


E  . p(f(j))p(g(j))  , 
a  <  j  <  b 


(14) 
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sm  b!' 


and  Sf  by 


S.  ,  -  E  e(hf(j)  +  £g(J))  , 

h  ,i  a  <  j  <  b 


S  _  *  £  min  Cl,  “"""""“rr^TTTTT}  3 

f  a  <  j  <  b  2«t||f(J)|| 


Lemma  4:  t  >  1 


1_  r  £  1  h<^l  +  — L.  £  £ 

2  1  <  h  <  t  1  <  /  <  t  2«2  1  <  h  <  c  1  <  £  <  C 


hi  -l  ,  ,5  c  j.  c 
hi  +  2  Sf  +  2  Sg  ' 


Proof:  Observing  tha t 


sin  2rthx  sin  2*2y  *  \  cos  2*(hx  -  iy)  -  \  cos  2it(hx  +  ly)  ,  (17) 


one  obtains  from  Lemma  3, 


p(x)p(y)  -  ^2  !  <  g  <  t  1  <  £  <  t 


cos  2x(hx  -  fa 
hi 


_1_  L  s  cos  2rt(hx  ±  ly 

2x2l^^^t  1  <  *  <  c 


I  "Uin  <*>  +  ■ln  (1’  mW1  • 
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Replacing  x  by  f(j),  y  by  g(j),  and  summing  with  respect  to  j,  Equation  (18) 
yields 


2*2  1  <  h  <  1  1  <  i  <  t  hi 

1  E  £  — 

2*2  1  <  h  <  t  1  <  2  <  t  hi 

18(D)  +  f  «,  ♦  i  V 


.  ?  .  c°8  2*(hf(j)  -  ig(j))  - 
a  <  j  <  b 


Z  cos  2x(hf ( j)  +  (19) 

a  <  j  <  b 


Thus 


Since 


<  •  y  £ 

2n2  1  <  h  £  1  1  <  i  <  t  hi 


>n2  1  <  h  <  t  1  <  i  < 


E  hi 


*g(j)) 


+  S,  +  S  . 
f  8 


cos  2it(hf (j)  -  ig(j)) 
a  <.  j  <  b  I 

.2  cos  2jt (hf ( j)  +  (20) 

a  <  j  <  b 


+ 


E  .  cos  2n(hf (j)  -  ig(j)) 
a  <  J  <  b 


a  <  J  <  bcos  2n(hf (J)  +  ig(j)) 


(21) 


the  lemma  follows. 
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) 


For  the  sequence  given  by 

Xj  -  j2),  (a,m)  -  1  , 

(22) 

is  necessary  to  estimate 

R  “  jl0lp4  j2  •  05  P<i^+T>2  '  P>  > 

(23) 

which  a,p  satisfy  0  <  a  <  1,  0  <  fJ  <  1.  Let 

f(j)  -  J  J2  -  a  , 

(24) 

g(j)  -  ^(j-*)2  -  P  ; 

(25) 

then, 

Lemma  5: 


h,i 


<  'J  2m(h+i,  m) 


Proof:  One  has 


h,  l 


Let 

d  -  (h,  £,m)  , 
h’  ■  h/d  , 
i'  ■  £/d  , 
m'  *  m/d  ; 


(26) 


\ 


(27) 
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f 


Chen 


h  ,£ 


i*?1  .ah'  .2  .  ai*  ,.,T%2N 

j£oe<~  i  +  —  <J+T>  ) 

m’-1  .ah'  .2  a V  ..  .2. 

jSo  *  +  ~  <*">  5 


Let 


then 


n'-l  ah»  2  at.'  2 

s'-jSo  ■ 


'tiftOi’  +  r))‘  +  JS=~:-  J) 

j-0  m'v  J  m'  J/ 


One  has 


,9  Bl’-l  m'-l  a  2  2  2a#<T 

Is’!2  -  k£o  +  l')(J2  -  k2)  +  ijf-  (j-k))  , 


and 


is-!2  -  S1  "&*•  4*** +  *’>“2  -  k2> +  ^  <J-k» 


Let 


j  «  k  +  v  , 


In  which  v  Is  a  new  summation  variable;  then 


s’l2  -  %l  +  i‘)vk)e4r(h.  +  i-)v2  +  2ifl  V) 


(28) 


(29) 


(30) 


(31) 


(32) 


(33) 


*  (34) 
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) 


and  hence, 


S*|2  <  K1^( h'+l')vk) 


(35) 


Let 


B  ■  (h •  +  i',m')>  b  *  “ — -  ^  •  ®"  m 


6  * 


B  ’ 


(36) 


then 


+  r>vk)  -  vk)  - 6  vk)  •  (37) 


m"-l 


Thus,  one  obtains 


i  i 2  m'_1 

IS* I  <8  VS0 


m  * - L  7ab 

&  rt> 


(38) 


By  direct  summation,  one  has 


kSQ1*^  Vk)  “  °’  m"  ^  2abV 


m",  m"  I  2abv 


Since 


(39) 


(ab,m")  -  1  , 


(40) 


m"  | 2abv  at  most  23  times,  and  hence, 


!  *  I  2  <  252m" 


25  m"  -  2Bm'  . 


(41) 
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Equations  (28),  (29),  and  (41)  now  yield 


h,/ 


<  d  •J  25m*  -  ^26d2m'  “  ^  26dm  -  s/  2m(h  +  i,m)  . 


(42) 


Lemma  5  yields  a  trivial  estimate  when  applied  to  S.  .  .  It  will  be 

n,  -n 

important  to  determine  an  accurate  estimate  for  this  quantity. 

Lemma  6:  1  <  t  <  >  v  <  h  <  t,  m  >  2t 


Sh,-h  -  0  • 


Proof:  One  has 


h,  -h 


J2-cj+->2» 


(43) 


Since  (a,m)  -  1,  the  sum  in  Equation  (43)  is  zero  when  m/2hT.  One  has,  for 
l<h<t,  »>2t. 


2  <  2ht  <  2tT  <  m  . 


(44) 


Thus,  m|2hT  and  consequently,  ■  0 


Lemma  7:  t  >  1,  m  >  2t,  1  <  t  < 

'Tm  _ 


R  < 


J2n2  1  <  *>  <  t  1  <  f  <  t 


^  V  (h  +  l,  m)  + 


V T1 


-Till  1  <  h  <  t  1  <  l  <  t 

h  4 1 


jjj  -  T,  m)  +  #  +f  S 


2  "f  T  2  g 
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Proof:  The  lemma  follows  immediately  from  Lemmas  4,  5,  and  6. 

Lemma  8:  For  t  >  3,  one  has 


E 

1  <  h  <  t 


3.27  In 


i<5<  tn  ■'<h  +  *>“> 

t  •  In  m) (2  +'/’2)v(m^  , 


(5.46  ln2t  + 
n 


in  which  v(m)  denotes  the  number  of  distinct  prime  divisors  of  m. 


Proof:  In  abbreviation,  define  S  as  follows. 


V m  j. 

■vT 2«2  1  <  h  <  t  1  <  £  <  t 


E  n/  ( h  +  2,m) 


(45) 


Let  h,£  be  restricted  so  that  (h  +  £,m)  ■  d;  then 


sT2x2  d|m  1  <  h  <  t  1  <  l  <  t  h* 


(46) 


(h  +  £,m)  ■  d 

The  set  of  integers  h ,1  for  which  (h  +  l, m)  ■  d  is  included  in  the  set  for 
which  d  |  h  +  l;  hence 


n/2  if  d|n 


S  <  ?  '/’d  E  si.. 

1  <  h  <  t  l<i<th^ 

d|h+i 


(47) 


All  integers  fall  into  d  residue  classes  modulo  d.  Let  r  denote  the  integers 
0, 1,2,. . .  ,d-lj  then  when  h  belongs  to  the  residue  class  represented  by  r,  i 
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belong*  to  the  class  in  which  d-r  is  a  member.  Let  Sr  denote  the  following 
sum; 


S  “  £ 

r  1  <  h  <  t 


h*r(mod  d) 
/*d-r(mod  d) 


(48) 


and  S'(d)  denote 


S’(d)  -  £  £  ~  ; 

1  <  h  <  t  1  <  £  <  t 


(49) 


then 


d  h  +  f 


<  s'<d)  , 

<F\ 2k  d  I m 


(50) 


»’«>  ’  S5  Sr  ' 


From  the  symmetry  of  the  sum  in  Equation  (48),  it  follows  that 


S  -  S  .  , 

r  d-r  ’ 


(51) 


and  hence, 


S'(d) 


S  +  2 
o 


1  <  rL<  d/2Sr  ' 


(52) 
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Setting  h  ■  cd  and  l  -  ed  in  which  c  >  1,  e  >  1  are  new  independent 
summation  variables,  one  has 

■(  z  ,  ^  '  ,N2 

d 


So  "  !2  (  1  <  cS<  t/d  c)  -  d2  (  1  <  c  <  t  c 


(53) 


Since,  for  t  >  3  , 


£  c<l+  f  -  ■  1  +  In  t  <  2  In  t. 
1  <  c  <  t  c  Jl  x  ' 


(54) 


one  obtains 


so<±JT1 

°  d2 


(55) 


For  r  satisfying  1  <  r  <  —  ,  let  h  ■  cd  +  r  and  £  *  ed  -  r;  then  h  >  1  implies 
c  <  0,  £  >  1  implies  e  >  1,  and 


Sr  —  o  <  c  <  cd+r  1  <  e  <  ed-r 


(56) 


One  has 


*  t-r  -i -<I+r 

<  c  <  __  cd+r  _  r  J0 


t-r 

d  dx  1,1.  t  .  1  In  t 

xd+r  r  d  r  -  r  d 


(57) 


Also,  one  has 


r  _J_  .  L  _ L_  <  i  £  _ L 

1  <  e  <  -pp  ®d+r  0  <  e  <  ~p  -1  ed+d"r  _d0<e<~S-le+| 


,  (58) 
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in  which  the  inequality  r  <  j  was  used;  further, 

0<.<V-1  ±iS2+l 


-  “  ~  d 


-J 


<  1.  3  +  In  t  <  2.  3  In 


Thus,  one  has 


r  1  <  2.3  In  t 

i  ^  ^  t+r  ed+r  ^  d 

i  <  e  <  — 


Equations  (56),  (57),  and  (60)  yield 


_  .  2.3  In  t  .  2.3  In  t 

Sr  <  rd  +  “ 


Equation  (52)  now  takes  the  form 


4  In2  t  .  2.3  In2  t  4.6  In  t  _  1 

d  d  l<r <§  t 


S>(d)  <  '  — ”  -c  + 

d* 


Since 


and 


4  In2  t  .  4  In2  t 
d2  £  d  ' 


in  which  the  inequalities  In  2  <  •  7  and  d  <  m  were  used,  one  obtains 


t-  (59) 


(60) 


(61) 


(62) 


(63) 


(64) 
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«> 


t  +  4.  6  In  t  *  In  ro 
'  *  d  d 


(65) 


Thus,  Equation  (50)  yields 


S  <  -~r  (7. 68  In2  t  +  4. 6  In  t  •  In  m)  Z  -pi- 

J2«2  dTm«7d 


(66) 


Let  0(a)  be  a  multiplicative  function  of  the  integral  variable  a,  and  let 


Q1  “2  °v 

a  "  P1  '  P2  Pv 


(67) 


be  the  canonical  factorization  of  a;  then  the  following  identity  holds  : 


J  0(d)  -  n  ti  +  e(p)  +  o(p2)  +  ...  +  ©(p0)]  . 
d|«  p|a 


(68) 


< 


Employing  Equation  (68)  with  0(d)  ■  1/ */' d,  a«*m,  one  obtains 


f  - r~  ■  n  (1  +  4-  +  ~—r — z  +  .. .  +  — f — -  ]  , 

d|m  v  d  p|m  P  (vp)^  ( vp)a 


(69) 


and  hence, 


Jf.75<p|n,'I+^+7^  + 


- 


pin  *P  -  1 


(70) 


Since  p  >  2,  one  has  vT p/(  >Tp  -  1)  <2+^2,  and  hence, 


df«775  (2+'r2)V<”>  • 


(71) 
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Equations  (66)  and  (71)  now  yield 


S  <  (7.68  In2  t  +  4.  6  In  t  •  In  m)(2  +  J~2)V^ 

'll* 


(72) 


Finally,  use  of  the  inequality  1/ n/~ 2  <  .  71  yields  the  result  of  the  lemma. 


Lemma  9:  For  t  >  3,  one  has 


•f n  j, 

<f 2n2  1  <  h  <  t  1  <  £  <  t  w 
h  +1 


Z  r~  '/(h  -  £,m)  <~  (8.52  In2  t  + 


2.84  In  t  •  In  m)(2 


Proof:  In  abbreviation,  define  S  by 


S  - 


•J~ m 


s[ 2  *2  1  <  h  <  t  1  <  £  <  t 

hj-t, 


L  -J  {h  -  £,m) 


(73) 


As  in  the  proof  of  the  preceding  lemma,  one  may  write 


S  “  ?  '/"d  2  2  r--  . 

/2r2  dTm  1  <  h  <  t  1  <  £  <  t 

(h  -  £,m)  *  d 

h^£ 


(74) 


Define  S’(d)  by 


S*(d)  -  2  2  7$  ; 

1  <  h  <  t  i  <  £  <  t 

”d|h  -  r 

h +1 


(75) 
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(  ) 


Chen,  since  Che  class  of  incegers  satisfying  (h  -  £,m)  -  d  is  included  in  Che 
class  d { h  -  i}  one  has 


™ 

When  h  belongs  Co  Che  residue  class  represenced  by  r,  Chen  i  belongs  Co  Che 
same  residue  class.  LeC  denote  the  following  sum: 


S  *  £  £  —  • 

r  1  <  h  <  t  l<2<th^' 
h  =  r  (mod  d) 

£  =  r  (mod  d) 
h  £ 


(77) 


(  » 


then 


S’(d) 


d£1 


’0Sr 


(78) 


Setting  h  ■  cd,  £  -  ed,  the  inequalities  1  <  h  <  t,  1  <  £  <  t  imply  1  <  c  <  t/d, 
1  <  e  <  t/d,  and 


S 


o 


L 

<  c  <  t/d 


In2  t 


(79) 


For  1  <  r  <  d,  let  h  -  cd  +  r,  l  ■  ed  +  r;  then  1  <  h  <  t,  1  <  £  <  t  imply 
0  <  c  <  0  <  e  <  and 


I  ) 
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Sr“o<c<iH  0<e<-^  (cd+r> <ed+r>  ‘ 
““  d  d 


(80) 


cf  e 


Due  to  the  symmetry  of  the  sum  in  Equation  (80), 


S  <  2  £  t_r  -~— 

r  —  o  <  c  <  c  jr  cd+r 


1  <  e  <  ■“  ed+r 
—  —  a 


(81) 


Equation  (57)  yields 


0  <  c  < 


— ! —  <  -  +  ln  t 
t-r  cd+r  -  r  d 


-  ~  d 


(82) 


further, 


E  “A?  <  4  £  7  <  4  (1  +  In  t)  <  2  ™  , 

1  <  e  <  Ijl  ed+r  d  1  <  e  <  t  e  d  d 


(83) 


and  hence. 


S  < 


4  In  t  4  In  t 


r  rd 


,2 


(84) 


Since  1  <  d  <  m  and 


d-1  1 

E.  —  <  1  +  In  m  , 
r-1  r  ' 


(85) 


Equations  (78)  and  (84)  yield 


S ' (d)  <  12-  -^n — -  +  -  .  1  ' — iil— EH 


d 


(86) 
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From  Equations  (76)  and  (86),  one  now  has 


S  <  (12  In2  t  +  4  In  t  •  In  m)  E  1 


/2s 


d  m 


/d 


(87) 


Equation  (71)  and  the  inequality  1/  /2  <  .  71  yield  the  result  of  the  lemma. 


Lemma  10:  (a,m)  *  1,  0  <  a  <  1,  0  <  p  <  1,  m  >  2nt 


c  0V(®)  m  In  m 

Sf  <  3. 2  “  , 

v(m)  m  In  m 


S  <3.2 
g 


Proof:  Consider,  in  the  sum 


m-1 


Sf  "  j£0  min  <1'  |(£  j2  _  a||>  ' 


(88) 


regrouping  the  terms  so  that  all  j  for  which 


a  j  ss  r  (mod  m) ,  0  <  r  <  m 


(89) 


constitute  the  sum  E{;  then 


m-1 

Sf  “  rS0  Sr  * 


(90) 


Let 


a  ai 


m  ■  2  p^  . . . .  Pk 


(91) 
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3 

be  the  canonical  factorization  of  m;  then,  by  a  standard  theorem,  the  number 

2 

of  j  for  which  a j  =  r  (mod  m)  is  T  ...  where  is  the  number  of  solutions 

2  aa  2  —  (X. 

of  aj  as  r  (mod  )  and  T  is  the  number  of  solutions  of  aj  =  r  (mod  2  }.  Also, 

one  has^ 

T£  <  2,  1  <  l  <  k  (92) 

T  <  4  . 


Hence, 


T  T,  ...  T.  <  4. 2 
1  k  — 


(93) 


and,  since  v(m)  ■  k  +  1, 


IT,  ...  T,  <2.2 
1  k  — 


v(m) 


(94) 


Thus, 


(95) 


and  hence, 


Sf  <  2.2V(m)  %  min  (1,  ) 


115  -  41 


(96) 


The  sum  in  Equation  (96)  will  be  estimated  by  consideration  of  four  cases. 

Case  1:  -1<—  -q<-t. 

-  m  -  2 


One  has 


||^  -  dj|  »  1  +  -  -  or  , 


(97) 
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and 


£  ,  min  (1,  - ” - )  <  1  + 

<  m(a-i)  2*t(l  +  -  -  a) 


0  <  r  <  m(a-|) 
m(a  -  -|) 


f 


min  (1 


2*t(l  +  -  -  a) 


-)dx, 


P L  \ 

<  1  +  m  J  min  (1,  j^j)du  * 


(98) 


<  1  + 


o 

m  In  itt 
nt 


Case  2: 


•i  <  -  -  a  <  0 

2  m  - 


One  has 


Hr  -  cm  =  a  - 1, 


(99) 


and 


met 


£  min  (1,  “ 

0  <  r  <  ma  2nt(a 


— <  1  +  / 
a  -  — )  ^  o 


min  (1,  - ”)d*  > 

o  2*t(a  -  -) 


1  +  m  f  min  (1,  j~)du  • 
o 


(100) 


<  1  + 


m  In  2rtt 

ut 


Case  3:  0  <  r  -  ct  <  • 

■-  —  m  <l 


One  has 


IS  '  ^  "  m  '  “  ' 


(101) 
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i  n  i 

„  E  „  min  (1,  - 7 - )  <  1  +  /  min  (1,  - - - )dx  , 

<  r  <  m  2*t(-  -  a)  ^ma  2«t(-  -  a) 

'm  ^  m 

<  1  +  m  f min  (1,  2ntU> du  ’  C102) 

j  o 


<  1  + 


m  In  2nt 


Case  4:  t  <  “  *  u  <  1  . 

— —  2  m 


One  has 


||-  -  a||  -  1  -  -  +  a  , 
Hm  m 


T.  min  (1,  - - ')  <  1  + 

i  ^  .  t  *  * 


m(a  +  — )<r<m 


2nt(l  -  -  +  a) 
m 


p  m 

^m(a  + 


min  (1, 


2nt(l  -  -  +  a) 

m 


-)dx  , 


pl  1 

<  1  +  m  j  min  (1,  7^)du, 
J  o 


<  1  + 


m  In  nt 


One  has 


S,  <  2.2v(m)[4  +  2m  IS-ill  + 

f  Kt  7TC 


and 
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Sp<8.2v(m)[l+^L^] 

r  rtt 


(106) 


Since  m  >  2irt,  one  has  In  2«t  <  In  m,  and  hence, 


Sf  <  8.2V<m)[l  +2Li£J2]  . 


(107) 


Also,  since  m  In  m  >  2nt,  one  has 


s  <  l2.2v(m)  EJijn 

f  rtt 


(108) 


Using  the  inequality  n  >  3  yields  the  result  of  the  lemma.  For  the  sum  S^,  one 
has 


m-1 

sg  “  j£0  min  (1> 


m-1 


- r - )  ■  L  min  (1,  - r - ) 

llr(J+T) 2-p)II  j  \\z  J  -ell 


(109) 


and  hence,  the  estimate  obtained  for  above  applies  also  to  S^. 
It  is  now  possible  to  state  the  first  main  theorem. 

THEOREM  1;  (a,m)  -  1,  m  >  36,  l<t<Vm,  0<a<l,  0<P<1, 


m-1 


j2-a)p(^(j+r)2.p)|  <\fm  [|  (2  +'/"2)v(m)  ln2m  +  30.2v(m)  In  m] 


,v(m) 


Proof:  Lemmas  7,  8,  9,  and  10  yield 


|R|  [13.98  In2  t  +  6. 11  In  t  •  In  m](2  +  «/2)v(m) 


+  15.2 


v(m)  m  In  m 


(HO) 


Since  *  >  9.85,  one  also  has 


0 
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| R |  c/m  [1.42  In2  t  +  .621  In  t  •  In  m](2  +  /2)V(m) 


+  15.  2V^  m  In  m  . 


Choose 


(111) 


t  -  •|\Tm  j 


(112) 


then,  since  In  t  <  j  In  «  , 


|R|  </m  (2  +sT2)V^>  In2  m  +  30.2v(m)  In  m] 


(113) 


The  conditions  t  >  3,  m  >  2*t  are  both  met  by  the  condition  m  >  36.  Since 
m/2t  -  ^m,  the  condition  1  <  t  <  m/2t  becomes  1  <  t  <'/m.  The  theorem  is  now 
established. 

The  autocorrelation  function  of  the  sequence  x.  »  {—  J2}  is  obtained 

J  ® 

immediately  from  Theorem  1  by  setting  a  ■  0,  p  ■  0,  and  recalling  that  i|r(T)  ■ 
R/m.  Hence,  one  has 

THEOREM  2:  (a,m)  ■  1,  m  >  36,  1  <  f  </m 

=0  |t(t)|  <  m2[|  (2  +^2)v(m)  In2  m  +  30.2v(m)  In  m]  . 


Consider  the  simultaneous  Diophantine  inequalities 
0  <  {f  J2)  <  Q,  0  <  (f(j-K)2)  <  P,  0  <  J  <  m  . 


014) 


Let  the  number  of  solutions  of  the  Inequalities  be  designated  by  T(a,f3),  then 
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0<cr,f»  ■  <ll5> 

a  2  fl  2 

ia  the  Joint  distribution  function  of  the  sequences  (-  J  },  l“(J+T)  )•  If  the 
sequences  (jjj  j2),  (^(J+t)2)  were  independently  equidlstrlbuted  over  (0,1),  one 
would  have  G(a,0)  ■  00.  It  is  the  present  object  to  determine  the  deviation  of 
G(a,0)  from  the  desired  Joint  distribution  00.  For  this  purpose,  let 

Ha(x)  -  a  +  p(x)  -  p(x  -  C0  ;  (116> 

then  Hq(x)  is  a  periodic  function  with  period  1  and,  within  the  initial  period. 


Hq(x)  -  1,  0  <  x  <  a  , 
•  o,  a  <  x  <  1  . 


(117) 


In  view  of  the  above  properties  of  HQ(x),  the  enumeration  T(a,0)  is  given  by 


(118) 


T(a,0)  “  ft  V^ft  • 

Lemma  11:  T(a,0)  -  00®  +  0  jSo^  ^  “  P  j§oP^m  ^  + 

a  %>&^ft  +  a  jSp<s<^>2  -p> + 

%<  ft  -  %<  J*  ■  a)  o&^ft  • 

%<  ft  P<Su*02  -  *>  +  J2  -  a>  -  p)  • 


Proof:  Use  Equations  (116)  and  (113). 
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_12:  (a,  m)  -  1,  m  >  36,  1  <  t  <Jm,  0  <  a  <  1,  0  <  P  <  1, 


m-l 


|T(a,P)  -  oe»l  <  2  I  jgp(J  J2)l  +  I  j£op(t  -  °9l  +  I  JZ-  p)l  + 


m-l 


a  .2 


m-l  a  2 

r  _  /**  i* 


'Jfm[|  (2  +  '/2)v(m)  In2  m  +  120  •  2v(m)  In  m]  . 


Proof:  Theorem  1  enables  the  estimation  of  the  sums  of  products  of  the  p- 

functlons  to  be  effected.  Also  observing  that 

J2>  >  <ll9> 

J2-  »  -  Sp^)2  “  »  a  (120) 

the  lemma  follows. 

In  order  to  estimate  the  sum  of  the  p-functions  in  Lemma  12,  the  following 
lemmas  are  required. 


Lemma  13:  (a,m)  »  1,  0  <  a  <  1, 


=*>  |“gje(^  J2  -  ha)  I  <  •J  2m(h,m)  . 


Proof:  Let 


(h,m) 

h' 

m* 


d  , 
h 

d  ' 


m 

d  * 


(121) 
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and 


s  -  J2  -  M  ; 


(122) 


then 


One  has 


M  ■  I$«<tJ2>'  -“I  Al«<^  Al 


,  ,o  7  m'-l  m'-l  hio  7  h'a  2 

l«l  A  ]£oe<^J  -^k> 


7  m'-l  m'-l+k  u.  2  h'a  2 

j2  A  j£k  ^  k> 


Introduce  a  new  summation  variable  v  by 


(123) 


(124) 


j  -  k+v  ; 


(125) 


then 


Isl2  -  <2  *  «£**>  S  -2  A  I  *«&*<»» 


m’-l  .  m’-l 


By  direct  summation,  one  obtains 


kfc1*(^ivk)  "  °»  m'/2h’av  » 

■  m',  m ' | 2h'av  . 


(127) 


Since  (h'a,m')  -  1,  m'|2h'av  if  and  only  if  m'|2v  which  may  occur  at  moat  twice. 
Hence, 
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I  S| 2  <  2d2m'  -  2md  . 


Tha  lemma  follows  on  taking  the  square  root. 


Lemma  14:  (a,m)  ■  1,  0  <  Q  <  1, 


(128) 


J2  -  “>l  <TS  !  <  &1J-.  b'r<M  +  ^  1»  » 

Proof:  Use  of  Lemma  1  yields 


.  'Tim 


-  -2 


(129) 


m-1 

+  r  min  (1, 

j®0  2  o_.na  ,2 


2nt||-  ?  -  a|| 


-)  • 


Setting  t  ■  ^'fta  and  using  Lemma  9,  one  obtains 

1  l"1;1., ha  .2. 


■m_1  ,a  ,2 

-  u;  |  v 

1  i  i  2 

The  lemma  now  follows  on  employing  Lemma  13. 


lj5oP<i  y  -  c£>\  <  E  1  r  £  Ij£d«(T  j2)l  +  6/m2V(m>ln  m.  (130) 

1  <  h  <  tt  v  m 


Lemma  15:  (a,m)  ■  1,  0  <  a  <  1,  m  >  36 


|^Jp(J  i2  -  a)l  <  3.92*/m  (2  +vT2)v(m)  In  m 


Proof: 


One  has 
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—  L  ^  5  l  r  £  - 

*  l<h<r>/nh  d|“  1  <  h  <  -'/m 


(h,m)  -  d 


E  sfd  Elr{ 

*  d|m  1  <  h  <  |N/m 

d)h 


Let  h  -  cdj  then 


•/im  _  r.  y  I  <  S’  r  JL  £  -  < 

*  d | m  l<h<T>/mh  *  d|m  ^  l<c<|'TmC 


d  |  h 


Use  of  Equation  (71)  now  yields 


*/m(.  15  +  .27  In  ■)  £  A 
d|m 


(131) 


(132) 


S’  £  £«J(h,m)  </m(.  15  +  .  27  In  m)(2  +  '/2)v(“) 


«  1  <  h  <  ■j  >/n  h 


Since  m  >  36,  one  has 


(133) 


.  15  +  •  27  In  m  <  .  312  In  m  , 


(134) 


and  hence, 


’S  z  (h,m)  <  .  312'/m  (2  +  */2)V^  In  n  .  (135) 

*  l<h<i/m 


Lenina  14  now  yields 
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J2  -  00  |  <  .  312  •>/m  (2  +*>/2)V^m^  In  m  +  6>Tm2v^  In  m  •  (136) 

Since  2/(2  +^2)  <  1  and  v(m)  >  1,  one  has 

.312  (2  +«/2)v(o)  +  6.2v(m)  <  3.92  (2  +'/2)v(n,)  .  (137) 

This  proves  the  lemma. 

THEOREM  3:  (a,m)  -  1,  m  >  36,  1  <  t  <-fm  , 

=0  |G(a,p)  -  op|  <  m2  (2  +sT2)v(m)  (|  In2  m  +  88  In  s)  . 

Proof:  Lemmas  12, and  15  provide  the  following  inequality. 

J. 

|G(a,e)  -  op  I  <  m  2[|  (2  +/2)v(m)  In2  m  +  (15.68(2  +'/2)V(m)  + 

(138) 

120.  2V(m))  In  m]  . 

Using  the  inequality 

2V(m)  <  .6(2  +'/2)V^m)  ,  (139) 

The  theorem  follows. 

a  2  fl  2 

Theorem  3  thus  demonstrates  that  the  sequences  (~  j4),  (“(J+t)*)  are 
approximately  independently  equidistributed  over  (0,1). 
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